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Abstract 

We investigate deuteron formation in nuclear matter at finite temperatures 
within a systematic quantum statistical approach. We consider formation 
through three-body collisions relevant already at rather moderate densities 
because of the strong correlations. The three-body in-medium reaction rates 
driven by the break-up cross section are calculated using exact three-body 
equations (Alt-Grassberger-Sandhas type) that have been suitably modified 
to consistently include the energy shift and the Pauli blocking. Important 
quantities are the lifetime of deuteron fluctuations and the chemical relaxation 
time. We find that the respective times differ substantially while using in- 
medium or isolated cross sections. We expect implications for the description 
of heavy ion collisions in particular for the formation of light charged particles 

at low to intermediate energies. 
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I. INTRODUCTION 



The formation of deuterons is the simplest multichannel reaction process in a heavy ion 
collision. Within a systematic quantum statistical approach the formation is driven by the 
Boltzmann collision integral that in turn requires the proper break-up cross sections. Besides 
electromagnetic disintegration/formation of the deuteron [np ^ 'jd) that is a generic two- 
body process, the treatment of three-body collision {NNN ^ Nd) is more involved, and 
requires a solution of the three-body dynamics in medium. Collision rates are a central input 
into modern microscopic simulations of heavy ion collisions such as the Boltzmann-Uehling- 
Uhlenbeck (BUU) or the quantum molecular dynamics approach (QMD) For the 
latest developments and recent applications to multifragmentation in heavy ion collisions 
see Ref. 0. Danielewicz and Bertsch started the numerical solution of coupled BUU 
equations including the deuteron formation. Later this model has been enlarged to include 
the production of three-particle clusters [0. 

So far these approaches are fed with experimental cross sections, e.g. Nd NNN for 
the deuteron break-up, and the medium effect is taken into account by the Pauli-blocking 
of final states. However, in principle, as shown earlier those cross sections are changed 
substantially, when medium effects are consistently included Besides including these 

medium dependent cross sections in the respective simulations that is currently on its way, 
it is an interesting and important question to analyze the change of time scales induced by 
medium dependent reaction cross sections. 

The correlated many-particle system is treated within the cluster mean field approach ||T^ 
that goes beyond the quasiparticle picture and in turn has been applied to a wide vari- 
ety of many-particle physics problems. This approach is intimately related to the Dyson 
Equation Approach to Many-Body Green functions |[TT|Jl^ . Both approaches lead to a sys- 
tematic decoupling of the equation hierarchy and thus to effective few-body equations. For 
the two-particle case these are known as Feynman-Galitskii or Bethe-Goldstone equations. 
Recently, we have derived a corresponding equation for the three-body system for finite 
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temperatures based on the formulation of Alt, Grassberger, and Sandhas (AGS) [|T3[. 

For zero temperatures three-body in-medium equations have been given in Ref. |jl2| for the 
related 2p-lh problem. 

The respective one-, two-, and three-body equations are derived in an independent par- 
ticle basis. That means we do not consider correlations in the surrounding matter that 
would eventually lead to a self-consistent treatment. However, such a treatment is techni- 
cally involved and only very recently progress has been achieved towards this direction by 



including two-hody correlations into the one-body spectral function ||15l. Since presently our 
aim is different, i.e. to extract the formation and equilibration time scales, and the densities 
considered are rather low compared to the ones that require a self-consistent treatment, and 
since many technical problems still need to be tackled we presently do not consider a such 
an approach for the three-body case. 

The life time and the relaxation time are determined by linearizing the Boltzmann equa- 
tion or the corresponding rate equations. Therefore the respective integrals can be evaluated 
in thermal and chemical equilibrium. To this end we assume only small fluctuations, i.e. 
linear response of the system. This will be explained in the following section. The necessary 
three-body reaction input on the basis of Faddeev type equations will be given in Sect. |T|. 
The derivation of these equations is based on the Green function approach, e.g. ||16|. In 



Sec. 13 we will present numerical results and in the last section we summarize and give our 
conclusion. 



II. REACTIONS 

The quantity of interest in the quantum statistical approach is the generalized quantum 
Boltzmann equation for the nucleon /Ar(p, r), deuteron /^(P, i?), etc. Wigner distributions. 
For the time being we assume symmetric nuclear matter. These distributions are driven by 
a coupled system of Boltzmann equations, see, e.g. Ref. 0, 

dtfN + d,U ■ drfN - drU ■ djj, = ^nU N , fd] (1 - f n) " ICTUn, fd] fw, 
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dtfd + dpU ■ dnf, - dnU ■ dpf, = V^[f^, /J (1 + /,) - JC^Vn, U U (1) 

where U denotes the total mean field potential. The coupling between the different species 
is through the collision integrals /C[/Ar,/d] that are merely shown for the deuteron case in 
the following two equations relevant in the present context, 



+ Jd^kJ d?kA \{kP\U\hk2)\%a^NafN{kut)fd{h,t)fr,{k,t) 
+ J d^k J d^kid^kid^ks \ {kP\UQ\kik2k^)\l^N^^^ 

X fNiki, t)/iv(A;2, t)fNik3, t)/iv(A:, t) 
+ ..., (2) 

where we used the abbreviations /tv = (1 — /at) and fd = (1 + fd), the ellipsis denote further 
possible contributions involving four and more body collisions (e.g. tp ^ dd, hp ^ dd) or 
processes like np ^ jd, and /CJ^'^*(P, t) is given by 
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/C7*(P,t) = Jd^kJ d^kid^k2d^k3 \{kik2k3\Uo\kP)\l 

X fN{ki, t)fN{k2, t)fN{k3, t)fN{k, t) 

+ ... (3) 

The quantity Uq appearing in Eqs. (0) and (H) is the in-medium break-up transition 
operator for the Nd NNN reaction and is calculated using the AGS equation as given in 
the next section. For the isolated three-body problem Uq determines the isolated break-up 



cross section c^[;^ via ||T3l|l7 



= ~ I d%d%d% \{kP\Uo\k,k2k3)\^ 

\vd -vn\S\ J 

X 2TT5iE' - E) (27r)35(3) + + k^) , (4) 

where \vd — vn\ denotes the relative velocity of the incoming nucleon and deuteron. 

So far the strategy has been to implement the experimental cross section into Eqs. (0) and 
This has then been solved, for a specific heavy ion collision 0. Using experimental cross 



sections respectively isolated cross sections may not be sufficient in particular in the lower 
energy regime. The cross section itself depends on the medium, e.g. blocking of internal 
lines or self-energy corrections of the respective three-body Green functions. This has been 
shown in Refs. [^|[ . While the inclusion of the medium dependent deuteron break-up cross 
sections in a coupled BUU simulation of a heavy ion collision is a congruous application, 
it is instructive to study the influence of medium dependent cross sections in a relaxation 
time approximation, in addition. This would enable us to already arrive at conclusions 
on the possible effect for typical conditions at heavy ion collisions. After linearizing the 
Boltzmann equation it is possible to define a break-up time for small fluctuations of the 
deuteron distributions. To do so, we consider only the term given in Eq. in the collision 
integral. For small fluctuations Sfd(t) = fdit) — from the equilibrium distribution 
linear response leads to 

§-6UP,t) = -i .6UP,t), (5) 

where the "life time" of deuteron fluctuations that depends on the deuteron momentum P 
and the nuclear density n^- has been introduced earlier 0, 

><I%{h)f%{h)f%{h)f%{k^) 2vr5(i? - Eo), (6) 

which can be related to the break-up cross section using Eq. (^. 

Chemical equilibration results from break-up and formation reactions, according to the 
terms that contain the transition operator Uq in Eq. and (j^). The change of the total 
deuteron density is given by the integration of the deuteron momentum in Eq. (|l]). For a 
homogeneous system we can neglect the drift terms and only the reaction terms will change 
the total density, therefore we obtain the following rate equation 

jUdit) = -a{t)nj^{t)nd{t) + /3(t)n^(t), (7) 

where the introduced rate coefficient for disintegration a is given by 



a{t)= J SPd^k j d^kid^k2d^ks \{kik2k3\Uo\kP)\ 



dN^pnN 



X Mk,,t)Mk2, t)Mks, t))f-^^l^^, (8) 

and an analogous expression holds for the formation rate f3. To obtain a time scale for the 
formation of chemical equilibrium, we consider a near-equilibrium situation. 

In chemical equilibrium both rate coefficients are constant and related by detailed balance 
(time reversal invariance) resulting in 



Here we denote the equilibrium distribution functions by n% and n% Assuming now small 
density fluctuations, viz. = n% + dn^ and Ud = n^^ — 671^/2, and after insertion into 
Eq. (|l]) a chemical relaxation time may be defined by 

1 5ndit) = ^ 5ndit). (10) 

at T,c\{nN) 

Using the definition of the rate coefficient a given in Eq. (P) leads to following expression 
for the chemical relaxation time 

r-i{n%) = f d'kd'P f%{k)fd{P)\vd-VM\ crUk,P) !^!k±^^ (n) 

The deuteron break-up cross section is defined as in Eq. (|^) but may now also include 
medium modifications of the transition operator Uq as will be explained in the next section. 
The deuteron equilibrium density nj] is related to the nucleon density via /x^ = 2yUAr for the 
chemical potentials. 



III. FINITE TEMPERATURE THREE-BODY EQUATIONS 

The basis to derive effective in-medium few-body Green functions for finite temper- 
atures is provided by the cluster mean field approximation. This approximation allows 
us to truncate the hierarchy of Green function equations in a controlled way. A general 



6 



introduction may be found in Ref |Ty]. This approach is related to the Dyson Equation Ap- 



proach [|TT],|T2[ extended here to finite temperatures |I^p!8| . It has been proven useful in the 



context of strongly coupled systems provided, e.g., in nuclear physics [jT2[. A calculable form 
is achieved by introducing ladder approximation. We consider generic elementary two-body 
interactions V2 only. The resulting equations for the decoupled one-, two- and three-body 
Green functions at finite temperatures (utilizing the Matsubara technique to treat finite 
temperatures [Q) will be given in the following. The one-particle Green function reads 



G^iz) = R?\z) = (z - e,)-\ (12) 
In mean field approximation the single quasiparticle energy ei is given by 

S^^(l) = Y: 12) - V2{12, 21)] /2, (13) 

2 

and /2 = 7^^(62) the Fermi function^ This is given by 

f(uj) = -—A (14) 

where f3 denotes the inverse temperature and n the chemical potential. The equation for 
the two-particle Green function G2{z) reads 

G2iz) = N2RP{z) + eS\z)N2V2G2{z), (15) 

Note that we use a matrix notation for the Fock space indices (1, 2, 3, . . .) for convenience. 
The matrix for the two-body resolvent R2'\z) is then given by 

R^^\l2,l'2'-z) = ^^^^ (16) 

Z — £i — £2 

and the Pauli blocking factor N2 by 



^Here and in the rest of the paper we use equilibrium distribution functions only, and therefore 
we use f% f for notational simplicity. 
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iV2(12,l'2') =(5ir522'(/l/2-/l/2) 

= 5n'522'(l-/i-/2). (17) 

We use the abbreviation / = 1 — /. The respective equation for the three-particle Green 
function relevant to describe three-body correlations in a medium is given by 

G,iz) = NsRf\z) + Rf\z)[N2V2]G,{z), (18) 

where we have introduced the notation 

3 

[N2V2]{123, 1'2'3') = Y.{N2V2Y^\123, 1'2'3'), (19) 

7=1 

{N2V2Y'Hl23, 1'2'3') = (1 - /i - /2)V^2(12, l'2')533'. (20) 

The last equation is given for 7 = 8. The channel notation is convenient to treat systems with 
more than two particles [|l^ . In the three-particle system usually the index of the spectator 



particle is used to characterize the channel. The Pauli factors A^3 and the resolvents R^^ 
read respectively 

iV3(123, 1'2'3') = 5iv522'hy{fif2h + /1/2/3) (21) 

i?f (123, 1'2'3'; z) = . (22) 

z — ei — 62 — S'i 

Note that [N2V2] ^ [iVal/z]^ and N^Rf = R^^N^. Although we assume a fermionic system 
the proper symmetrization of the equations is treated separately. 

If the correlated pair, e.g. '12', and the spectator particle, e.g. '3', are uncorrelated in the 
channel (3) we may define a channel Green function G^^\z). We generalize to the channel 
(7). The channel Green function is then defined by 

Gp\z) = z^E ^G2{u;,) G,{z-u;,). (23) 

The summation is done over the bosonic Matsubara frequencies ux, A even, ux = ■n\/ {—i[3) + 
2fi. The equation for the channel Green function is derived in the same way as for the total 
three-particle Green function given in Eqs. ([I8|). The result is (no summation of 7) 
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G^^\z) = NsRf\z) + rI'\z){N,V,Y^^G^^\z). (24) 
We are now ready to define the proper transition operator Uai3, 

Gs{z) = 5^pGi''\z) + Git\z)U^p{z)Gf\z), (25) 
which in the zero density hmit coincides with the usual definition of the transition operator 



with the correct reduction formula to calculate cross sections [T^. After proper three-body 



algebra we arrive at the following equation for the transition operator in medium, viz. 

N,U^p{z) = (1 - 5„^)[i?f (^)]-i + Y.iN2V,Y''^G^\z)U,^iz). (26) 

For iVa 7^ the equation may be multiplied by N^'^ from the left. We now connect this 
equation to the two-body t matrix. To this end we define the following transition channel 
operator T^"'\z), 

G^^\z) = N,rI>\z) + Rf\z)Ni''^Tt\z)N,RPiz). (27) 

The definition has been chosen so that is leads to the standard Bethe-Goldstone equation in 
the two-nucleon subchannel, T3 (z) = T2{z)633', where z is the frequency of the two-body 
subsystem embedded in the three-body system. Inserting eq. ( pTf ) into eq. (|^) leads to the 
subchannel Bethe-Goldstone equation 

T^^\z) = Vi''^ + {V2N2Y'^'>RP{z)Ti''\z). (28) 

With the help of this equation it is possible to connect the AGS type equation directly to 
the two-body sub t matrix 

NsU^^iz) = (1 - 5^p)[R^^\z)]-' + Y: Nt^^Ti^'\z)Rf\z)N,U,p{z). (29) 

Here we have written the Pauli factors occurring due to the surrounding matter explicitly. 
Note, that through eq. (HD T^^^ {z) is also medium dependent. 

1/2 1/2 

Since > Q me may multiply Eq. (^) by N-^ form the left and by A'"3 from the 
right, and introduce t/*^ = Nl^'^U^pNl^'^ . Then 
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= (1 - (^)]-^ + E N^'^'Ni^^Tt\z)Nl/'Rl>\z)U;,{z). (30) 

We may now introduce the transition channel operator T-^^'^'' = N^^^"^ N2'''T^'^'^ nI^"^ so that 
the equation looks formally as for the isolated case. The respective channel t matrix equation 
then is 

The explicit form of the effective potential arising in this equation reads (e.g. in the (3) 
channel) 

1/3*^^^(123, 1'2'3') = A/'3"^/^(123)A/'f ^(123)\/i^^(123, l'2'3')A/'3^^(l'2'3') 

= (1-/1-/2)^/^(1-/3 + ^7(51 + 52))-^/' 

x\/2(12, 1'2')533.(1 - /3 + g{ev + e2')f'\l - fy - hf'^ (32) 
^ (1 - /i - /2)'/V2(12, 1'2')(1 - h> - f2'f'\ (33) 

The last equation holds for P <^ f only and g{uj) = [e^^"^"^^) — 1]^^. Utilizing this ap- 
proximation eq. (^) has been solved numerically using a separable potential for the strong 
nucleon-nucleon interaction ^j. 



IV. RESULTS 

Since the Green functions have been evaluated in an independent particle basis the one-, 
two-, and three-particle Green functions are decoupled in hierarchy, as given in Eqs. (|12|), 
(p!5|), ([T8|). To solve the in-medium problem up to three-particle clusters the one-, two- 



and three-particle problems are consistently solved. This leads to the single particle self- 
energy shift Eq. (pr3D, the two-body input including the proper Pauli blocking Eq. (|3TD, and 
eventually to the three-body scattering state. 

For technical reasons we have used some reasonable approximations valid at smaller 
effective densities. The nucleon self-energy has been calculated via Eq. (|TB|), but in the 
three-body code we have approximated the self-energy by effective masses, i.e. 
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For simplicity and speed of the calculation, we use angle averaged Fermi functions {{N2)) 
and ((iVg)), e.g. 

{{N3)) = j^-^ J dcos9g d cos 6p d(j)qd(j)p N3{p,q, P^.m.), (35) 
where the angles are taken with respect to Pc.m. and p, q are the standard Jacobi coordi- 
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We have solved the three-body equation neglecting terms of the order p, using eq. (|33D 
as the driving kernel in the AGS equations, which is sufficient, since the calculation is done 
below the Mott density of the deuteron. 

A rank-one Yamaguchi potential has been used in the calculations of the deuteron 
break-up cross section that includes the coupled ^Si-'^Di and the ^Sq channels only. The 
parameters are given in Table |. To get an impression of the quality of the calculation the 
isolated cross section is given in Fig. |l| along with the experimental data on neutron deuteron 



scattering pl |. 



Unlike the isolated three-body problem the integral equation depends on the relative 
momentum of the three-nucleon system with respect to the medium, since the Fermi func- 
tions depend explicitly on the center of mass momentum of the three-body system. As a 
consequence one has to solve the three-body problem at a finite center of mass momentum 
-Pc.m. = ki + k2 + in a medium that may be considered at rest Pmod = 0. However, techni- 
cally it is more convenient to let the three-nucleon system rest Pc.m. = and the surrounding 
medium move with Pmcd = —(^1 + k2 + k^). This procedure results in the least change of 
the three-body algebra and is possible since the dependence on Pc.m. — Pmed is only through 
the Pauli blocking factors thus parametric. 

For small momenta Pc.m. and small relative energies of the dA^-system the infiuence of the 
medium is most pronounced. In Fig. ^ we see that the in-medium cross section is significantly 
enhanced compared to the isolated on. So it is not a justifiable approximation to include 
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medium effects only by the Pauli-blo eking factors of the final states (see eq. (0)). That would 
effectively decrease the break-up cross section and might lead to wrong conclusions. The 
dominant effect is the weakening of deuteron binding energy (Mott-effect). This provokes 
also a shift of the threshold to smaller energies with raising densities. We observe that for 
higher energies the medium dependence of the cross section becomes much weaker, which a 
posteriori justifies the use of isolated cross sections (along with the impulse approximation) 
when higher energies are considered 0. The dependence of the cross section on the center- 
of-mass momentum Pc.m. is shown in Fig. |^ for a fixed density and temperature. With the 
definition of a medium dependent cross section according to Eq. (^ we can rewrite for the 
deuteron break-up time Eq. (|^) 

^b"u (^) ='^Jdk^ \Vd - VNWy,^{k, P)f{kM). (36) 

To get an impression how much the medium dependence of the cross section shortens the 
deuteron break-up time we show in Fig. ^ the result of the above Eq. (|36| ) with use of the 
medium cross section a^^ calculated by solving Eq. (^) compared to the isolated a^^. For 
small deuteron momenta the break-up time is shorter by a factor of three at n = 0.007 fm^^ 
and T = 10 MeV.. 

We also calculate the chemical relaxation time given in Eq. (|Tl]). In Fig. ^ we see again 
that the medium dependence of the cross section shortens the time scale especially at higher 
densities. 

V. CONCLUSION AND OUTLOOK 

Our results show that medium dependent cross sections in the respective collision in- 
tegrals lead to shorter reaction time scales. Chemical processes become faster. This also 
effects the elastic rates that are related to thermal equilibration. 

The basis of this result is the cluster Hartree-Fock approach that in our approximation 
includes correlations up to three particles in a consistent way. The equations driving the 
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correlations are rigorous. The respective one-, two- and three-body equations are solved, 
in particular for the three-particle case Faddeev/AGS type equations have been given in 
Ref. [p|,p|,p!4| ■ The AGS approach is particularly appealing since it allows generalizations to 
n-particle equations in a straight forward way. Results for the three-body bound state in 



medium will be published elsewhere As expected from the deuteron case, the triton 
binding energy changes with increasing density up to the Mott density, where Et = 0. 

The production rates, distribution of fragments, spectra etc. of light charged particles in 
heavy ion collisions at intermediate energies may change because of the much smaller time 
scales induced through the medium dependence compared to the use of free cross sections 
(respectively experimental cross sections). 
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FIG. 1. A comparison of the total, elastic, and break-up cross sections nd — > nd, nd 
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with the experimental data of Ref. [21|. 
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FIG. 2. In-medium break-up cross section at T = 10 MeV and 
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= 0. The isolated cross 

section is shown as a solid line, other lines show different nuclear densities as given in the legend. 
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FIG. 4. Deuteron break-up time at T = 10 MeV and nuclear density n = 0.007 fm^^. Solid 
line using the medium dependent cross section as given in Fig. ^, and dashed line using the isolated 
cross section. 
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FIG. 5. Relaxation time for small fluctuations of the deuteron density from chemical equi- 
librium at a temperature of T = 10 MeV as a function of the nucleon density. Line coding as in 
Fig. I 
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TABLES 

TABLE L Parameters for the separable potential used to solve the three-body AGS equation. 
The form factors are guiip) = Ct^ip^ /{p^ + P'^iY^'^^^- The binding energy of the isolated deuteron 



is Ed = 2.225 MeV and of the isolated triton Et = 7.937 MeV H 



Co C*2 (3o j3i 

[MeV fm-^] [fm-i] 



^Si - ^Di 8.63237 -38.8017 1.24128 1.94773 
^So 8.87026 1.17328 
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